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Abstract: The magnetically suspended Control Moment Gyroscope (CMG) has the 
advantages of long-life, micro-vibration and being non- lubricating, and is the ideal actuator 
for agile maneuver satellite attitude control. However, the stability of the rotor in magnetic 
bearing and the precision of the output torque of a magnetically suspended CMG are 
affected by the rapid maneuvers of satellites. In this paper, a dynamic model of the agile 
satellite including a magnetically suspended single gimbal control moment gyroscope is 
built and the equivalent disturbance torque effected on the rotor is obtained. The 
feedforward compensation control method is used to depress the disturbance on the rotor. 
Simulation results are given to show that the rotor displacement is obviously reduced. 

Keywords: Control Moment Gyroscope; dynamic modeling; feedforward compensation; 
magnetic bearing 



1. Introduction 

With the development of high resolution imaging on Earth, the requirement of satellite agility is 
increasing. The magnetically suspended Single Gimbal Control Moment Gyroscope (SGCMG) is 
chosen as the ideal actuator for its large output torque to perform rapid maneuver tasks by using the 
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magnetic bearing instead of traditional mechanical bearing [1-3]. It has the advantage of high 
precision, long-life, micro-vibration and being non-lubricating [4]. It can be applied for agile maneuver 
satellite attitude control. For agile satellites, the slew rate is often at the level of l~10°/s [5], but when 
the satellite slews fast or the gimbal rotates rapidly, the rotor displacement will be increased and the 
precision of output torque will be decreased, even leading to rotor instability, so research on agile 
satellite attitude control using magnetically suspended SGCMGs is very important. 

For the control of magnetic bearing systems, scholars have reported some research results [4-7]. 
For example, cross feedback control of magnetic bearing systems [6], filtered-X least mean square 
(FXLMS) algorithm for moving-gimbal effects [4], and the compound control for moving-gimbal 
effect compensation to increase the response of the gimbal system [7]. The research is all based on a 
static base, and the influence of rapid satellite maneuvers on the rotor displacement is not analyzed. 

For the satellite attitude control based on a Control Moment Gyroscope (CMG), there are some 
research results [8-14]. Reference [12] gives a model of the CMG mechanism and the CMG output 
torque in body frame. In Reference [13], aiming at solving the problem of high precision and high 
stability attitude control of satellites, the influence of magnetically suspended rotor dynamic imbalance 
and static imbalance on the vibration of the base, and the stability of a magnetically suspended control 
system are investigated. The disturbance torque on the base resulting from the imbalance vibration is 
reduced. In Reference [14], a dynamic model of a spacecraft with a magnetically suspended single 
gimbal control moment gyroscope is given, but the slew rate of the satellite is limited at 10 rad/s level 
(approximate 0.57s). How to use magnetically suspended SGCMGs for the agile maneuver attitude 
control of the satellite is not mentioned. 

In this paper, the dynamic model of the agile satellite including a magnetically suspended single 
gimbal control moment gyroscope is built and the equivalent disturbance torque effected on the rotor is 
obtained. The feedforward compensation control method is used to depress the disturbance on the rotor. 
Simulation results are given to show that the rotor displacement is obviously decreased. 

2. The Satellite Attitude Dynamic Model Including a Magnetically Suspended SGCMGs 

Different from the mechanical CMG, a magnetically suspended rotor is used in a magnetically 
suspended CMG. The translational position of the rotor relative to the magnetic bearing can be 
controlled actively. Besides the high-speed rotation degree, the magnetically suspended CMG has five 
more degrees of freedom than the mechanical CMG, and thus the modeling process is complex. 
Firstly, the coordinate frames are needed to be defined. In the following, n = 1, 2, 3, 4 denotes the n-th 
CMG. Figure 1 gives the relationship between inertial frame system, orbit frame, and satellite 
body frame. 

Fi: Inertia frame O^Y^. The origin is the Earth's core, points to the vernal equinox, while 
O.Z { points to the North Pole. 

F 0 : Orbit frame OXYZ . The origin is the centre of mass of the satellite, n x and n y represent 

O 0 0 0 *~ 7 J 0 0 O O A 

the roll axis and the pitch axis, respectively, both locating in the orbital plane. 

Ft,: Satellite body frame O b X b Y b Z b . The origin is the centre of mass of satellite. o b x b , o b Y b 

and o b Z b are along the principal axes of inertia of the satellite. 
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Figure 1. The relationship between the coordinate frames. 




Fc mg , n • Th e n -th CMG frame O cmg n X cmg n Y cmg n Z cmg n . Fixed with the satellite, but it is determined 
by configuration. O cmgn is the rotor centre of mass of the n -th magnetically suspended SGCMG, 
X n is the gimbal axis, Y is the spin direction of the rotor when the gimbal angle is in 

zero position. 

F g n : The n -th gimbal frame O g n X g n Y g n Z g n . Fixed with the gimbal, the frame can rotate with the 

gimbal. When the gimbal angle is in zero position, the gimbal frame coincides with the CMG frame. 
F f n : The n -th magnetic bearing installed frame O f n X f n Y f n Z f n . It is fixed with the gimbal frame. 

This frame is obtained by 45° rotation of the gimbal frame and about the Y g n axis. 

F r n : The n -th rotor frame O r n X r n Y r n Z r n . It is fixed with the rotor. This frame does not spin with 
the rotor. O f n coincides with O r n when the rotor is not suspended. 

Cp' : denotes the coordinate transformation matrix from frame ¥ x to frame F 2 . 

In Section 2.1, the rotor dynamic model of single gimbal magnetically suspended CMG is analyzed. 
For convenience, the index n is omitted. For example, F g represents the gimbal frame O g X g Y g Z g • 

2.1. The Dynamic Model of Magnetically Suspended Rotor 

The Euler equation is used to build the magnetically suspended rotor dynamic model. A magnetically 
suspended SGCMG consists of the rotor and the gimbal. The magnetic force is produced from the 
unique relation between the rotor and the bearing. In order to perform the satellite attitude control, the 
direction of angular momentum is changed through magnetic bearing torque by the gimbal rotation, 
and then the gyroscope torque is generated, and is transmitted to the satellite by the magnetic bearings. 

In this paper, the rotor is supposed to be asymmetric, and the dynamic and static vibration induced 
by the geometry shape are ignored. Only the force induced from the gimbal rotation of the CMG frame 
relative to the inertial frame is considered. The rotation relationship between CMG frame and the 
gimbal frame is given in Figure 2. 
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Figure 2. A sketch of the magnetically suspended SGCMG. 




As shown in Figure 2, 8 , 6 and ^ are the angular velocity of CMG frame relative to gimbal frame. 
8 is the rotation angle from CMG frame to gimbal frame. 8 is the gimbal rate about X c mg axis. 
Figure 3 denotes the view from A-side about the Y cm ^-direction in Figure 2. It shows the relationship 
between the gimbal frame and the installed magnetic bearing frame, fx and f z denote the radial 
magnetic bearing force about the X-axis and Z/-axis. f Y is the axial magnetic bearing force, a and /? 
are the rotation angle of the rotor about the X/-axis and Z/-axis, respectively, a and (3 are the angular 

velocity about the X/ -axis and Z/ -axis, respectively. Q is the rotor speed relative to the rotor frame, 
and it is about Yy -axis, ijx and i AZ are the magnetic bearing control current about X/ -axis and 

Z/-axis of the A-side. 



Figure 3. Sketch map of the magnetic bearing installation. 



A-side ^ X 



45» "T^X, 




Firstly, only one magnetically suspended SGCMG is analyzed. M r is the torque acting on the rotor. 
By using the Euler equation, the magnetically suspended rotor dynamic model in the rotor frame can 
be obtained: 



M r =H r +<xH r 



(1) 
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In Equation (1), (o r ir denotes the absolute angular velocity in rotor frame. It includes <o r/ , o icmg 9 
and <o g . <o£ is the angular velocity of the rotor. It is the magnetic bearing frame relative to the rotor 
frame, and is shown in the magnetic bearing frame. (o^f g is the angle velocity of CMG frame relative 
to the inertial frame. <o g is the gimbal angular velocity: 

o r =C f ((o f f +C f C 8 (o cmg ) n\ 

ir f \ rf g cmg ig ) \A) 

The absolute angle velocity Cl r ir of the rotor includes the spin rate H r and the rate (o r ir of the 
rotor frame: 

o; = or + < = or + c> « + c;c*x« ) (3) 

Because a and /? are very small, then cosa ~ 1, cos /? ~ 1, sina ~ a, sin/? a/? ~ 0. Then: 



C r = 



cos j8 sin /? 0 
-sin/? cos/3 0 
0 0 1 



1 0 0 

0 cos or sin a 

0 -sin or coscir 

cos 45° 0 -sin 45 

0 1 0 

sin 45° 0 cos45 c 





~ 1 


p 


0" 




or 




-p 


1 


a 




0 




0 


-a 


1 




A 



(4) 



By rotating 8 about X cmg axis from the CMG frame to the gimbal frame, then: 





~1 


0 


0 




~d 




'f 




~8~ 




S 


c g = 

cmg 


0 


cosS 


sin<J 


, CO* = 


0 


cmg = cmg g = 
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-sin S 
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0 







(5) 



where § , 6 and (j> show the absolute angular velocity of the CMG frame, and the detailed expression 

will be given in Section 2.2. By substituting the coordinate transformation matrix and the relative 
angular velocity, then: 



^-|<j-^cos£+#sin c>j + + cos c> + ^sin 



Q -—/?(<£- 0cos£+#sin d\- d/3 + (^0cosd +^smd^-—aid + ^cosS -dsinS) + 



^-(s+<f>cosS-0sin £>] + /?-or(#cos£+^sin 



2 \ 

In Equation (1), H r is the rotor angular momentum in the rotor frame: 

"/„ 0 0 



(6) 



H =f or 



0 Iry 0 

o o r 



(7) 



where l r denotes the wheel moment of inertia in the spin direction. I rx and I ry are the radial inertia of 
the magnetically suspended rotor in the x and y direction, respectively. I rz is the axial inertia of the 
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rotor. The angular momentum variation rate relative to the rotor frame can be obtained from the 
time derivatives: 

- ^ 'cos + <f)8 sin + 0 sin 8 + 68 cos 8\ 



H r = 



+j30 cos £ + /?# cos £ - pOS sin 8 + p(/)sinS + pip sin £ + /fyk? cos S + a 

-pS-pS + p^cosS + p^cosS-p^hinS-pOsinS-pOsinS-pdScosS 
+6c8 +a 'S + d^cos 8 + or0 * cos 8- a0S sin 8 -dO sin S-aO sin 8 -aOS cos 8 
+£l + 6cos8- 68 sin £ + <p > sin £ + (j>8 cos 8-ap + ap 

+ (/) cos S - $8 sinS - 0 sinS- OS cos 8^ 

dO cos S-aO cos 8 + sin 8- d(j)s\n 8- aip sinS-a&S cos 8+ 0 



(8) 



By substituting Equations (4)-(8) into Equation (1), the rotor torque in rotor frame can be obtained: 



M r = 



^11 i 



2 \ 

V 



|i+^cos J-^sin£j + /? +(/ n -7 ry )(^cos^+^sin^) ^-|i+^cos J-#sin£j + /? 



c) - cos £ + ^sin 8 + # sin 8 + ftj cos J J + fid cos 8 + fi(/)sin8 + a 



(I rx -I rz ) f<? + 0cos£-0sin + £-^cos£ + #sin£ +rir 

V 2 A 2 



— ^- > 5^ + ^cos^- > 5^sin^ + ^+^cos^-dr^sin^ 



V 2 



+Q + # cos 8 -08 sinS + <p sin S + (f)ScosS-ocfi + afi 



2 \ 



</>cos8+6s'm8 \ + cc 



+ (^ry -/ K )((9cos J+^sin J) — |<j- ^cos <9sin + a 



2 V 



<p cos 8-08sin8-6 sin 8-68 cos J - cos 8-CC(j)sm.8+ (3 



(9) 



From Equation (9), it can be seen that, for every magnetically suspended SGCMG, I^QS is the gyro 
coupled torque produced by the gimbal rotation relative to the inertial frame. IS is the inertial 
coupled torque produced by the gimbal. I cos 8 and I ^£10 sin S are the gyro coupled torque 
produced by the satellite rotation. I n ip cos 8 and I n 0sinS are the inertial coupled torque produced by 
the satellite rotation. I^P and I^Qd are the gyro coupled torque produced by the rotor tilt relative 
to the magnetic bearing. IJx and I n fi are the inertial coupled torque produced by the rotor 
accelerating tilt relative to the magnetic bearing. 

Equation (9) includes the gimbal angular velocity and the satellite angle acceleration. It shows that 
there is strong coupling between the rotor, the gimbal and the satellite. Moreover, there are many 
trigonometric functions, for example, sin 8 and cos 8 , which shows strongly nonlinear characteristics. 
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2.2. The Torque of Magnetically Suspended SGCMG Clusters 

In this paper, the pyramid configuration using four magnetically suspended SGCMGs is adopted 
(Figure 4). The frame of the w-th CMG can be obtained by rotating the satellite body frame y n about 

the Z b axis, and then rotating <J n about the axis, y = [y 7 2 7 3 n] T =[ 90 ° 180 ° 270 ° °] T , 

G = [a, a 2 a 3 <J 4 f =[53.13° 53.13° 53.13° 53.13°f. 

Figure 4. The pyramid configuration installation of four magnetically suspended SGCMGs. 




By rotating the rotor frame 45° in the negative direction of the Y g axis, rotating -8 about the X g 
axis, and then transforming according to the installation matrix of the pyramid configuration, the 
attitude transformation matrix C b r from rotor frame to satellite body frame can be obtained. Because 

the magnetic bearing gap is small, the rotation displacement can be ignored relative to the rotation of 
the gimbal and the satellite, then C f r ~ E : 



r cmg g f r cmg g f 



where 



cL = 



sin ex sin 7 cos y -cos asm y 
-sin a cos y siny cos a cos y 
cos a 0 sin a 



r^cmg _ ( f^g V 1 _ 



cos 45° 0 sin 45° 

0 1 0 
-sin 45° 0 cos45 c 



The transformation matrix from gimbal frame to body frame is: 

sin a sin y cos ycos 8 - cos a sin ysmS - cos ysinS- cos a sin ycos 8 
- sin a cos y sin ycos 8 + cos a cos ^sin 8 - sin ^sin 8 - cos a cos ^cos 8 
cos G sin a sin 8 sincrcos^ 



r-\b s~ib r-\cmg 

g ~ cmg g ~ 



(10) 



(11) 



(12) 
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Equation (10) can be written as: 



c = 



sfl V2 sfl V2 V2 V2 

— sin £ cos y + — sin o sin 7 + — cos £ cos <r sin 7 cos £ cos 7 - sin 8 cos cr sin y — — sin 8 cos y + — sin a sin 7 - — cos £ cos cr sin y 

V2 . . V2 . V2 V2 . . V2 . 72 

— sin 0 sin 7 sin (J cos 7 cos <) cos cr cos y cos 0 sin y+ sin 0 cos cr cos 7 sin 0 sin 7 sin cr cos 7H cos o cos cr cos 7 

2 ' 2 2 / / ' 2 ' 2 / 2 / 



V2 V2 . 

cos cr cos 0 sin cr 

2 2 



sin £ sin cr 

Then, the absolute angular velocity of the CMG frame can be obtained 



V2 V2 . . 

— cos cr + cos^>sin<j 

2 2 



(13) 





sin (J sin y 


— sin cr cos y 


coscr 






~8 




cos^ 


sin y 
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at 

wy 
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-coscr sin y 


cos a cos y 


sin cr 






0 



(14) 



The torque generated by the rotor in body frame can be obtained after the sum of four CMGs: 



(15) 



By Equation (15), M 6 , the torque of the magnetically suspended SGCMG cluster acting on the 
three inertial principal axis of the satellite can be obtained. 

2.3. Satellite Attitude Dynamic Model Including Magnetically Suspended SGCMGs 

The satellite attitude dynamic equation including magnetically suspended SGCMGs can be obtained 
by the law of angular momentum conservation: 

J©* + < x ( J©* ) + M b = u d (1 6) 

where J is moment of inertia of the satellite. (o b ib is the angular velocity vector relative to the inertia 

frame, and it is shown in the body frame, is the environmental disturbance torque, such as 

aerodynamic drag gravity gradient, solar radiation pressure and Earth magnetic torque. <°^ X ( J(0 ^) 

represents the gyroscope torque generated by the satellite rotation. 

Equation (16) can be simplified to J(h b ib =u c +u d , where u c is the command control torque acting 

on the satellite. The command torque for magnetically suspended SGCMGs is: 

Ug = -u c - ©* xH g -<x(j<) (17) 

where H g is the total angular momentum of the actuator relative to the satellite body frame. The 
gimbal angle of four magnetically suspended SGCMGs is § = S 2 S 3 S 4 ] T , where S n represents the 
gimbal angle of the /2-th CMG. By substituting y into Equation (12), the total momentum of four 
magnetically suspended SGCMGs is: 

- cos a x sin 8 X - cos S 2 + cos cr 3 sin S 3 + cos S 4 



cos 8 X - cos cr 2 sin S 2 - cos S 3 + cos cr 4 sin S 4 
sin a x sin 8 X + sin cr 2 sin S 2 + sin cr 3 sin S 3 + sin cr 4 sin S 4 



(18) 



where K - I ry Q denotes the rotor angular momentum. Then the variation rate of CMG angular 
momentum in the satellite body frame can be obtained: 
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u,=H,=C(8)8 



(19) 



where C(8) is: 



C(8) = 



-cosc^ cos<5j 

-sin^ 
sin o x cos S x 



sin£ 2 
-cos<7 2 cos S 2 
sin <j 2 cos 8 2 



cos <J 3 cos £ 3 



-sin^ 



sin S 3 cos <J 4 cos S 4 
sin cr 3 cos £ 3 sin <J A cos £ 4 



(20) 



The command gimbal rate is calculated by the robust pseudo-inverse steering law with null 
motion [15]: 



8 = 8 C +8 AT 



(21) 



where 5 C denotes the command value calculated by robust pseudo-inverse steering law, 5^ denotes 

the command value calculated by null motion. 

The attitude kinematics are described with the quaternion of the satellite body frame relative to the 
orbit frame, and then the attitude kinematics equation is: 



• 1 ^ 
q = -q®co o , 



(22) 



where q = [q 0 q x q 2 q 3 f represents the attitude quaternion, <g> means quaternion multiply. The 
satellite body rate relative to the orbit is (o ob = & b - Ag> 0 . A is the attitude rotation matrix relative to the 
orbit frame. o o = [0 -a> 0 of is the orbit velocity. 

3. Feedforward Compensation for the Rotor Disturbance Torque 

When the satellite maneuvers quickly, the satellite angular velocity and gimbal rate will introduce 
extra torque on the rotor, and the rotor tilt motion will be enhanced. The disturbance torque can be 
compensated by the feedforward control. Figure 5 gives a diagram of the feedforward compensation. 

Figure 5. Sketch of the feedforward compensation. 
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In order to keep the rotor displacement close to zero, the controller provides the control voltage to 
maintain the rotor stability. However, the rotation of satellite and gimbal introduces disturbances to the 
rotor. Feedforward compensation is used to reduce the influence on the rotor. The compensation does 
not change the control system stability. 

The disturbance caused by satellite and gimbal is given in Equation (9). When the satellite body and 
the gimbal are fixed, 8 , ^ and 0 are all zero. The torque is generated by the rotor. Equation (9) is 

written as: 



M r (0): 



I^a + lJ 



(23) 



Assuming that the rotor acceleration is zero, Q = 0 , the rotational inertia of the non-rotation axis is 



uniform, I G 



1^. By ignoring the first order item that includes a and /?, the above 



equation can be simplified as: 



M'(0) = 



0 



(24) 



By deleting M r (0) from M r , the equivalent disturbance torque acting on the magnetic 
bearing control system can be obtained: 



M r = 

LJM ~ext 



V2" 



72" 



I^Q^S + <pcosS-0sinSj + I rr j c ^ - 0 c< > 8 + (j)8 s i n S + 6 s in < > -!- 9( ) c< c ; 



0 



^^I^Q^8-(pcos8+0sm8^ + ^^I rr ^8+^cos8-(p8sin8-0sm8-08cos8^ 



4~2 



(25) 



Since 7 rz Q » I rn the gimbal and satellite have more influence on the rotor than their angle 
acceleration. It is to say, the gyro coupling torque is more obvious than the inertia coupling torque. In 
designing the steering law, the limited gimbal angular velocity must be considered to prevent the 
instability. The gimbal angle acceleration should be limited to improve the dynamic response ability. 



4. Simulation Results and Analysis 

The initial satellite attitude angle is [-20° 50° 30°f . The target satellite attitude angle is [o° o° o°f . 
The moment of inertia of the satellite is [12 12 6] T kg m 2 . The initial gimbal angle is 
[90° -90° 90° -90°f . The radial inertia of the magnetically suspended rotor is 
4 =/ rv = 4=00034kg m 2 9 the axial inertia is 4 = 0.0052kg- m 2 9 the speed of the rotor is 15,000 r/min. 
Q = 500;rrad/s . 
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4.1. The Comparison of with and without Feedforward Compensation not Using an Attitude 
Control Loop 

In this section, aiming at the dynamic characteristics of a single gimbal magnetically suspended 
CMG, the simulation is performed without the attitude control loop. The gimbal acceleration is 
120°/s 2 , the maximal gimbal rate is 10°/s, and the simulation period is 3 s. The gimbal rate limit is 
4iax =10%. Figures 6 and 7 are the simulation results of with and without feedforward compensation. 

From the local zoom in of Figures 6(a) and 7(a), it can been seen that, the vibration displacement range 
is decreased to some degree after using the feedforward compensation, and then a faster response and 
the higher precision of the output torque can be obtained. From the comparison of Figure 6(b) with 
Figure 7(b), it can be seen that the overshooting of control current is reduced by using feedforward 
compensation, and the power consumption is decreased. Because the magnetic bearing control 
frequency is highly relative to the satellite motion, the disturbance torque induced from the satellite 
angle velocity can be compensated by the feedforward loop. 

Figure 6. The curves without feedforward compensation, (a) Magnetic bearing torque; 
(b) Control current; (c) The rotor tilt angle; (d) The rotor tilt angle velocity; (e) The rotor 
displacement of A-side; (f) The rotor displacement of B-side. 
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Figure 6. Cont. 
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Figure 7. The curves with feedforward compensation, (a) Magnetic bearing torque; 
(b) Control current; (c) The rotor tilt angle (d) The rotor tilt angle velocity; (e) The rotor 
displacement of A-side (f) The rotor displacement of B-side; (e) The rotor displacement of 
A-side; (f) The rotor displacement of B-side. 
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Figures 6(c) and 7(c) give the extent of tilt of the rotor relative to the magnetic bearing. It can be 
seen that after using feedforward compensation, the tilt extent is reduced, but the relative tilt angle 
velocity is increased, which can be seen from Figures 6(d) and 7(d). Figures 6(e,f) and 7(e,f) 
correspond to the displacement of the two sides of the rotor. It can be seen that the rotor displacement 
in the maneuvering process is reduced apparently by the feedforward compensation. 

4.2. The Comparison of with and without Feedforward Compensation Using Attitude Control Loop 

The attitude control period is 0.5 s. The maximal gimbal rate is ^ max -10°/s. The simulation 
results of the three-axis attitude control of the satellite with feedforward compensation are given in 
Figures 8-13. Figure 8 is the satellite attitude angle. Figure 9 is the satellite angle velocity. Figure 10 
gives the gimbal rates of four magnetically suspended SGCMGs. Figure 1 1 gives the gimbal angles of 
four magnetically suspended SGCMGs. Figure 12 is the satellite angle acceleration. 

Figure 8. Satellite attitude angle. 
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Figure 9. Satellite angle velocity. 




Figure 10. The gimbal rates of magnetically suspended SGCMGs. 
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Figure 11. The gimbal angle of magnetically suspended SGCMGs. 
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Figure 12. Satellite angle acceleration. 
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Figure 13. Singularity measurement curves. 
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The magnetic bearing force in each channel depends on the magnetic bearing coefficient, the 
magnetic gap of the two sides of the magnetic bearing, the bias current, the rotor displacement and the 
control current, etc. When the magnetic bearing design is finished, it is determined by the magnetic 
bearing coefficient and magnetic bearing gap. By adjusting the control current, the magnetic bearing 
force is changed to assure the convergence of rotor displacement, but the magnetic bearing force has a 
nonlinear relationship with the rotor displacement and the control current. The magnetic bearing 
controller is used to make sure the stability under the working point. 

In Figures 14-19, the left column gives the rotor X-Y displacement, the magnitude of rotor 
displacement and control current without feedforward compensation. The right column gives the 
results with feedforward compensation. 
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Figure 14. The rotor displacement without feedforward compensation. 
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Figure 15. The rotor displacement with feedforward compensation. 
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Figure 16. The magnitude of the rotor displacement without feedforward compensation. 
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Figure 17. The magnitude of the rotor displacement with feedforward compensation. 
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Figure 18. The magnetically suspended rotor control current without feedforward compensation. 
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Figure 19. The magnetically suspended rotor control current with feedforward compensation. 
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It can be seen that, after using feedforward compensation, the rotor displacement is reduced to 
about 30%, which decreases the risk of the rotor instability. Figures 18-19 show that the control 
current does not vary apparently. Although the rotor displacement and the given control current are 
reduced, extra control current is needed for feedforward compensation. It is to say, the command 
torque is invariant. The torque is output by the control current to assure the stable suspending of the 
rotor, and so the control current did not vary apparently. By using feedforward compensation, the 
maximal displacement of the rotor is less than 10% of the magnetic bearing protecting gap, which can 
ensure the requirement of the stability of magnetically suspended rotor. Figures 20-23 give oc ? (3 \ a 
and (5 of the magnetically suspended SGCMG rotor, which describe the motion of the rotor relative to 
the magnetic bearing. 

Figure 20. The magnetically suspended rotor tilt curves relative to rotor center before 
compensation. 
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Figure 21. The magnetically suspended rotor tilt curves relative to rotor center after 
compensation. 
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Figure 22. The magnetically suspended rotor angle velocity relative to rotor center before 
compensation. 
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Figure 23. The magnetically suspended rotor angle velocity relative to rotor center after 
compensation. 
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By using feedforward compensation for the magnetically suspended control system, the response of 
angle acceleration and angle velocity becomes faster owing to the increase of the control current, 
which decreases the rotor angle displacement. There is a magnetic saturation problem in the magnetic 
bearing design, so there exists a magnetic bearing control current saturation limit. When the satellite 
and the gimbal rate exceed the limit, the additional disturbance torque on the rotor is large, which 
exceeds the extent of the control current, and then the rotor becomes instable, even resulting in CMG 
failure. This will influence the stability of the whole system. When the maximal gimbal rate limit 
is^max = 12% , the results of not using feedforward compensation are shown in Figures 24 and 25. 
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Figure 24. The gimbal rate of four magnetically suspended SGCMGs. 
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Figure 25. The rotor X-Y displacement without feedforward compensation. 
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The terminative condition of the simulation is whether the rotor displacement exceeds the magnetic 
gap (100 |im) or not. The above figures show that, the simulation finish in less than 1 s. The reason is 
that magnetically suspended SGCMG3 becomes instable. Because the gimbal angular acceleration 
command limit is not performed, the overshooting of gimbal servo system results in the real gimbal 
rate being larger than 12% , and then the equivalent additive magnetic bearing torque is large, which 
results in an increase of the rotor displacement. There exists the magnetic bearing control current 
saturation limit, it is difficult to provide sufficient magnetic bearing torque, which results in the rotor 
displacement exceeding the approximately linear zone and reaching the magnetic gap instantaneously, 
and then the magnetically suspended rotor becomes instable. In the end, the stability of magnetic 
bearing control system and satellite control system are destroyed. 

The curves of gimbal rate and rotor displacement with feedforward compensation are given in 
Figures 26 and 27. 
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Figure 26. The gimbal rate of four magnetically suspended SGCMGs with feedforward 
compensation. 
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Figure 27. The rotor displacement with feedforward compensation. 




It can be seen that the rotor displacement is reduced by the feedforward compensation. Although 
there is overshooting of the gimbal rate, even exceeding 14°/s, the displacement of the rotor is less 
than 5 |im. This ensures the normal use of the magnetically suspended SGCMG. There are also some 
research results that are related to the compensation and disturbance attenuation [16-19]. In the future, 
the advanced algorithm will be employed. 

4.3. Analysis on the Output Torque 

The output torque of magnetically suspended SGCMG includes the inertial coupled torque and the 
gyro coupled torque, which is produced from the coupling of the satellite, the gimbal and the rotor. 
Traditionally, the inertial coupled torque is ignored when Equation (19) is used to calculate the output 
torque according to the Jacobian matrix and the gimbal rate. In this paper, the total torque is produced 
from the magnetically suspended rotor dynamic equation, as shown in Equation (15). By using this 
method, a more actual output torque of the actuator can be obtained. This is important for the analysis 



Sensors 2012, 12 



16985 



of the agile satellite attitude control system. Equation (15) also includes the strong coupling 
relationship between the rotor, gimbal, and satellite. This results in strongly nonlinear in the 
trigonometric function item that is relative with the gimbal angle. Figure 28 shows the torque 
calculated from Equation (19). Figure 29 gives the total output torque calculated from Equation (15). 

Figure 28. Torque calculated using Equation (19). 
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Figure 29. Torque calculated using Equation (15). 
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From the simulation results, it can be seen that the trend of torque curves calculated using 
Equations (15) and (19) are consistent. Figure 29 shows the influence of the rotor imbalance vibration 
and gimbal system, which is high frequency compared with the satellite attitude control system, so in 
the simulation, the torque model in Equation (15) should be used, which can emulate the dynamic 
characteristics of a magnetically suspended SGCMG more really. In the feedforward compensation for 
a magnetically suspended rotor, the compensation should be performed according to the rotation rate 
of the CMG frame with respect to the inertial frame, so in a situation of satellite rapid maneuver, the 
dynamic modeling of the magnetically suspended SGCMG is very important. The simulation results 
show that it coincides with the theoretical analysis results. 
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5. Conclusions 

With magnetically suspended SGCMGs mounted as a pyramid, the modeling of satellite attitude 
dynamic is built, and feedforward compensation control is used. The simulation results are given to 
show that by using the feedforward control, the displacement of the rotor is reduced. This is important 
for the future semi-physics experiment, and becomes the basis for the application of a single gimbal 
magnetically suspended CMG in agile satellites. 
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